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A non-equilibrium thermodynamics model able to analyze the combined effect of 
diffusion and adsorption in porous materials is proposed. The model considers the 
coupled dynamics of the diffusive phase, described by a diffusion type equation, and 
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the adsorbed phase governed by a generalized Langmuir equation. It is shown that 

o i 

the combination of diffusion and adsorption can be treated as diffusion in an effec- 
tive medium, with an effective diffusion coefficient depending on the characteristics 

X" 

of the adsorption kinetics. The interplay between these effects leads to the appear- 
ance of different regimes in which the effective diffusion coefficient exhibits peculiar 
behaviours such as a decrease due to the presence of the adsorption process and a 
non-monotonous behavior resulting from particle interactions. We outline applica- 
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tions of the model developed to diffusion in porous materials such as zeolites and 
bentonites, and to engineered highly porous materials. 
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1. INTRODUCTION 

Diffusion is a fundamental irreversible process intervening in the evolution 
of many out of equilibrium systems. At sufficiently large scales it is successfully 
described by Fick's law obtained from non-equilibrium thermodynamics^ and, 
at shorter time scales, by means of kinetic equations of the Fokker-Planck type 
derived from mesoscopic non-equilibrium thermodynamics^-. 

When the particles diffuse through a heterogeneous or complex medium, 
the structure and the physicochemical properties of the medium may give rise 
to the appearance of additional kinetic processes that affect the dynamics of 
the particles^— . An important example of this situation, with applications 
ranging from biomedicine to engineering, is diffusion in the presence of ad- 
sorption in which particles moving through the pores of the material have 
certain probability to become attached at the walls of the pores^i~— . Two 
possibilities can be analyzed: total or partial immobilization of the particles 
at the surface^ 1 ^ 1 ^ 1 ^. The kinetics is thus dominated by a set of transport 
equations that may present different relaxation time scales depending on the 
physicochemical properties of the surface and of the bulk fluicU^ 1 ^ 1 ^. 

Our purpose in this article is to propose a simple model able to analyze the 
mutual influence of diffusion and adsorption processes in micro and nanopores, 
a situation observed in zeolites, bentonites or in new engineered highly porous 
materials^ 1 ^^— . We are interested in provide a simple theoretical model able 
to explain, at least in qualitative form, how this mutual influence is responsible 
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for interesting effects observed in experiments^ and recently studied via molec- 
ular dynamics and cellular automaton techniques on the diffusion properties of 
the particles in the pore^— . Examples of these effects are the non-monotonic 
behavior of the diffusivity as a function of particle concentration or particle 
loading of the pore, or even the strong suppression of self- diffusivity as the 
number of particles in the pore increases^— . The model describes diffusion 
of one component in another at rest by means of the Fick law and adsorp- 
tion through a generalized Langmuir equation obtained from non-equilibrium 
thermodynamics^. Using the fact that the concentration of the diffusing par- 
ticles depends on that of the adsorbed particles, due to mass conservation, 
we formulate an effective medium theory in which the overall process is de- 
scribed as effective diffusion, with a diffusion coefficient depending on the 
characteristics of the adsorption process. Our analysis on the effects of ad- 
sorption on bulk diffusion is complementary to the usual theoretical approach 
to the problem which makes emphasis on the surface dynamics of the adsorbed 
phased— >±tt— .Here, we are considering the case of isothermal adsorption. The 
generalization to the case when non-isothermal effects, such as pervaporation 
associated to the mass transport through a zeolite membrane could also be 
performed by following a non-equilibrium thermodynamics analysis similar to 
that of Ref. [30]. 

The effective diffusion coefficient is calculated for a particular adsorption 
kinetics. It is shown that it includes a factor accounting for its deviation for the 
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case of absence of adsorption. At high concentrations of adsorbed particles, a 
virial expansion is proposed to estimate the effects of interactions. We conclude 
that the effective transport properties of the system strongly depend on the 
adsorption-desorption rate constants as well as on temperature through the 
corresponding virial coefficients. The great multiplicity of behaviours found 
reveals the strong influence that one of these proceses exerts on the other. 

The paper is organized as follows. In Section 2, we present the model and 
analyze the case of total immobilization of the adsorbed phase. Section 3 is 
devoted to study the effect of surface diffusion of the adsorbed phase on the 
diffusion/ adsorption process. Finally, our conclusions are presented in Section 
4. 

2. INTERPLAY BETWEEN ADSORPTION AND DIFFUSION 

Non-equilibrium thermodynamics provides a general framework to describe 
both adsorption and diffusion process and to analyze their mutual influence^. 

In many instances, diffusion and adsorption occur simultaneously. Particles 
in the bulk phase diffuse and arrive at the surface where they become attached 
due to the presence of molecular forces. This last process corresponds to 
adsorption and has been modeled by considering it as a diffusion through a 
potential barrier-. 

When the pore is very narrow, the particles practically move along the pore, 
in a thin central layer (see Fig. 1). Under these conditions transport is quasi- 
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FIG. 1. Diffusing particles in a quasi-onedimensional micropore. The orange circles 
represent the mobile particles in the porous at time t. The light grey circles repre- 
sent the mobile particles at time t that become adsorbed at time t + t. The lines 
schematically represent the trajectories they follow. The black circles represent the 
particles adsorbed at time t + r 

onedimensional and the corresponding descripton can be performed in terms 
of averaged concentrations depending on the coordinate along the axis of the 
pore. It is convenient to stress here that we will only consider the case when 
the diffusing particles have a coherent structure, such that their geometry can 
be characterized by an effective radius which is much larger than the molecules 
of the solvent. This assumption avoids to consider the rotational dynamics of 



the particles with more complicated geometries. 

After averaging the three dimensional diffusion equations over (y, z)- 
directions, it can be shown that the evolution equation for the average number 
concentration c\(x) in the bulk is given by^ 

where J\ is the diffusion current of particles having dimensions of number of 
particles per unit area and time along the ^-direction, and j has dimensions of 
number of particles divided by volume and time and represents the net amount 
of particles leaving the bulk by crossing the dotted lines indicated in Figure 
1. In similar way, the evolution equation for the concentration c 2 (x) of the 
adsorbed phase is 

-jH = -£*(*)-:>+*.. P) 

Here ^(x) is the corresponding diffusion current density and we used the 
fact that the net amount of particles arriving at the surface is jin — ~j and 
that jcs is the reaction term associated to chemisorption. In the following, for 
simplicity we will only consider the case of physical adsorption, that is, when 
j cs = 0. Further details on the derivation of Eq. (j2J) are given in Ref. [4]. 

The expression for the currents entering Eqs. and follows from the 
entropy production rate per unit mass a of the system^. This quantity is given 
by a = o"i + cr 2 , where the contribution associated to the bulk is 

1 T 9 f*\ 
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and the entropy production at the surface is given by^ 

where fii(x,t) (i = 1,2) represent the corresponding non-equilibrium chemi- 
cal potentials. It is convenient to mention that, by hypothesis, no chemical 
reactions occur in the bulk and therefore Eq. Q has no contribution due 
to chemical reactions. From Eqs. ([3]) and (J4]) and assuming linear relation- 
ships between fluxes and forces, one finds: J\ oc (dfii/dx), J 2 oc (cty^/cte), 
j OC jl\ /i 2 . 

In the dilute case at constant temperature, the last relation implies that 
the flow of particles from the bulk to the surface j becomes proportional to 
the pressure difference which in turn is proportional to the concentration (as- 
suming Henry's law in the bulk). This flow has two contributions, j = j a + jd, 
one for the particles leaving the bulk, j a , and one for the particles arriving 
to the bulk from the surface, jd- The first contribution is proportional to a 
constant rate k a giving the probability per unit time that the particles in the 
bulk is adsorbed by the surface. However, this probability must depend on the 
available surface. In Langmuir approach, this factor was assumed as propor- 
tional to the concentration of particles at the surface: 1 — c 2 /c™ ax with c™ ax 
the maximum of adsorbing sites at the surface. However, if interactions exist 
among these particles, then it depends upon the non-linear function ^(£2) 
which takes into account excluded volume interactions of the adsorbed phase 
at the surface^ . That is, k a = k a Z 2 (c2) with Z 2 = 1 — C2/c 2 nax for ideal systems 



and Z 2 (c 2 ) arbitrary in the general case. Hence, the flow j a can be written 
as: i) j a = —k a Ci(l — c 2 /c™ aa: ) for Langmuir and ii) j a = —k a CiZ 2 (c2) for the 
non-ideal case. In both expressions k a represents the adsorption rate constant. 
The second contribution jd is proportional to the covered surface and is usu- 
ally expressed as jd = —k d c 2 , where k d is the desorption rate constant. Using 
these relations in Eqs. ([T]) and one obtains 



d_ d_ 

dt ° l dx 



ox 



+ 3, (5) 



for the mass balance in the bulk, whereas for the adsorbed particles in the 
nonlinear case we have 

dc 2 



d d 
—c 2 = k a ciZ(c 2 ) - k d c 2 + — 



D 2 — 

ox 



(6) 



In writing Eqs. (5) and (6) we are considering the concentrations c\ and c 2 
in terms of volume. This allows us to use the mass conservation condition 
c = Ci + c 2 , in order to considerably simplify the description with out loss 
of generality, and takes into account the fact that surface diffusion is strictly 
a three dimensional process that occurs in a small volume surrounding the 
surface of the solid. Thus, Eqs. (EJ) and together with mass conservation 
c = ci + c 2 , describe the dynamics of the particles in the pore. It is conve- 
nient to mention that the adsorption rate constant follows an Arrhenius law 
in which this quantity depends on temperature Ref. 4. As a consequence, the 
interplay between adsorption and diffusion could be controlled by changing 
the temperature of the bath. 
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Instead of proceeding with the numerical solution of Eqs. (j5J) and (JH]), 
we will use mass conservation in order to adopt a simpler effective medium 
description for the interplay between adsorption and diffusion processes in 
which the boundary term j in Eq. (jSJ) does not appears explicitly. This can 
be done because this term is coupled to the reaction term in d^J) , which in fact 
determines a relation between the concentrations C\ and Ci corresponding to 
the adsorption isotherm. For instance, in the stationary case with no surface 
diffusion effects, j is strictly zero. However, even in this case, we now that 
the diffusion of particles in the bulk is modified by the adsorption process, in 
such a way that D\ = -Di(c 1; c^). Using the mass conservation relation and the 
corresponding adsorption isotherm, this effect can be incorporated through an 
effective diffusion coefficient Z^ e //(c).->^ 1 ^ The effective diffusion current for 
the total number concentration is 

Oc 

Note that the effective diffusion coefficient D e ff(c) can be associated to an 
effective activity a e ff(c) = cf e ff(c) in such a way that we may write 

D cff (c) = BTB [l + d -MM 

oln\c\ 

where / e //(c) is the effective activity coefficient and B is the mobility which 
in the present case can be well approximated by the Stokes friction coefficient, 
B = 6irr]R e ff where t] is the dynamic viscosity of the solvent and R e ff is 
the effective radius of the particle. Using equations ((Tj) and ([8]) in the mass 
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conservation law for the bulk yields 



dc d dc 
dt dx eff dx 



(9) 



It remains now to determine the expression of D e ff(c) which, after using the 
adsorption kinetics and the mass conservation condition, can be expressed as 
a function of the total concentration c. This will be performed by analyzing 
different cases in the next section. 



The case in which adsorbed particles are completely immobilized follows 
from the general model by imposing that the effective diffusion current J e ff, 
referring initially only to the mobile fraction c\ is given by 



The diffusion process can also be described by means of the effective Fick's law 
(0). This fact implies a relationship between the diffusion coefficients D\[c\) 
and D e fj(c) that can be established by considering the constraint c = c\ + C2, 
with the subindex 2 corresponding to the immobile fraction. Since c 2 depends 
on Ci through the adsorption isotherm: c 2 = 0(ci), this gives c = C\ + 0(ci) 
or equivalently c\ = Ci(c). Thus, using the equality dci/dx = (dci/dc) dc/dx, 
we finally obtain 



2.1. Total immobilization of the adsorbed phase 



£ ) i(ci,c 2 ) 



<9ci 
dx 



(10) 



D. // ( C )=A( C )D 1 (c) ) 



(11) 
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where the factor A(c) = dci/dc accounts for deviations of the diffusion coeffi- 
cient with respect to its value in the case of absence of adsorption. 

The influence of adsorption on diffusion enters through the evolution in time 
of the concentration of the immobilized particles c^. An interesting analysis 
on how distinct kinetic regimes may appear in adsorption-diffusion interplay 
is discussed in Ref.— . This dependence can be determined by means of the 
generalized Langmuir adsorption equation^ 



When adsorption occurs in time scales much shorted than the characteristic 
time of diffusion, T a & << Tp, we may assume that the adsorbed phase is 
practically in equilibrium. Under this condition, one has from (fl2~|) 



where a = kd/k a . For the case of a non-equilibrium stationary state in which 
dc-i/dt = Jads, with J ai is a constant current of particles, one obtains 



Low coverage: When the concentration of the adsorbed particles is small 
we may assume Z 2 — 1 — £02, with £ = l/c™ ax the inverse value of the total 
concentration of adsorbing sites c™ ax in the system^^. Then, Eq. (IT2*|) yields 
the Langmuir isotherm^ 




(12) 




(14) 
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Using Eq. (fT51) in the condition c = C\ + <f>(ci), one obtains 

ci(c) = ^-(c-a+y/(c-a- l) 2 + Aca - l) , (16) 

where c = c£ is the total concentration scaled by the total concentration of 
adsorbing sites. From Eq. ffl6|) . the correction factor A(c) becomes 

A(c) = lfi + V (it) 

2 y ^/(c - a - l) 2 + 4ca y 
In Fig El the correction factor A(c) is represented as a function of the total 
concentration c for different values of a in the equilibrium case. From now on, 
we will assume for simplicity that £ = 1cm 3 . Values of £ different from the 
unity only rescale the results, the qualitative behavior of the functions remains 
the same. When adsorption is dominant, as occurs for a ~ 0.01, the effective 
diffusion coefficient becomes negligible for concentrations c < 0.4cm -3 . In- 
creasing the importance of desorption, a ~ 1.0, makes A(c) more significant 
taking values of about 50% to 70% of Z?i(c). As expected, this tendency be- 
comes even more manifest when desorption is faster than adsorption, a » 1. 



In similar form, for the stationary case we have 



C2 = Ma) = - ~rr-—T — ■ (18) 

a + Ci rc a a + c x 



and therefore the correction factor A s (c) is in this case 



... 1/ c — a — 1 + 2a 

A a (c) = - 1 + , . (19) 



2 



c-a-l) 2 + 4(J + ca) 
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FIG. 2. The factor A(c) from (|17p in terms of concentration for different values 
of the parameter a. Solid lines correspond to equilibrium adsorption (|15p whereas 
dashed and dotted lines correspond to the non-equilibrium stationary case (|19p for 
J = 0.1cm -3 and J = 1.1cm -3 , respectively. The values of a used were (from red 
to purple) a = k^/ka = 0.01, 1.0, 3.0 and 30.0 and £ = 1cm 3 . 

here we have introduced the parameter J = J a ds/k a , having dimensions of 
number concentration. 

Dashed and dotted lines in Fig [2] show the stationary non-equilibrium val- 
ues of A(c) as a function of the total concentration c for the same values of 
a in the equilibrium case and for J = 0.1cm -3 (dashed) and J = 1.1cm -3 
(dotted). The stationary adsorption process is still dominant for a ~ 0.01, 
but in this case the effective diffusion coefficient takes values considerably 
larger than those of the equilibrium case, 10% for J a d s = 0.1cm -3 and 30% for 
Jads — 1.1cm -3 . As in the equilibrium case, this effect decreases when increas- 
ing a. This behavior of the correction factor A(c) as a function of particle 
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concentration in the bulk is similar to the one observed for the effective dif- 
fusion coefficient in experiments and simulations^ 1 ^ 1 ^. The relation between 
both quantities can be established by means of Eq. ( TTTT) . Unfortunately, a 
quantitative comparison is not possible due to the fact that reported results 
in simulations are given in terms of the so-called loading parameter measuring 
the number of particles in the pore. 

High coverage: When the concentration of adsorbed particles is higher, 
one has to take into account nonlinear contributions to Z(c 2 ) through the virial 
expansion 



where b is a constant having dimensions of square volume and may take pos- 
itive or negative values depending on the nature of the interactions between 
adsorbed- adsorbed particles and may also depend on temperature. It is con- 
venient to mention that the expansion (1201) is consistent with a second order 
virial expansion in which interactions (excluded volume and direct ones) be- 
tween adsorbed particles are taken into account. This is the case since the 
factor Z(c2) plays the role of a fugacity. More important is that Eq. f[2"0"j) 
allows us to go beyond the classical Langmuir model. Using this relation, we 
obtain the following adsorption isotherm 




(20) 




(21) 
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FIG. 3. The nonideal correction factor A n j(c) in terms of the total concentration c, 
for b = (1 + a) 2 /4 and different values of a = kd/k a = 0.01, 1.0, 3.0 and 30.0 (from 
red to purple) and £ = 1cm 3 . i) The solid lines represent the behavior of A n j(c) for 
equilibrium adsorption whereas dashed lines correspond to the ideal case. Interac- 
tions between the adsorbed particles introduce a nonmonotonic behavior of the effec- 
tive diffusion coefficient for a ^ 1, which depends on the value of b. ii) The nonideal 
correction factor for a = 2, £ = 1cm 3 and 46/(l+a) 2 = 0.2, 0.4, 0.6, 0.8 and 1.0cm 6 
(from red to orange). 

where the subindex ni indicates the non ideal case under consideration. Note 
that Eq. (1211) is valid for b > 0. Substitution of this equation into the mass 
balance leads to the relation 



c = C\ 



a + ci£ 
2bc 1 



l±Jl-4bc 2 1 /(a + c 1 2 



(22) 



This equation is a third degree polynomial equation for c\{c) that has a real 
solution from which the correction factor A„j(c) can be calculated. The re- 
sulting relations can be used to estimate the corrections due to the non-ideal 
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nature of the adsorbed phase. The condition b < (l+a)/2, obtained from Eqs. 
(J2TJ) and 022]), guarantees the existence of real values for the concentration c. 

Figure O) shows the nonideal correction factor A n j(c) in terms of the total 
concentration c by assuming £ = 1cm 3 . For a < 2, interactions reduce the 
magnitude of the effective diffusion coefficient although it increases monotically 
with c. However, for a ^ 1, interactions between particles of the adsorbed 
phase may cause a non-monotonic behavior of the effective diffusion coefficient 
for certain values of b, making it to decrease for larger values of c as it is shown 
in detail in Figure Elii) for which a = 2. This behavior implies that an optimal 
value for the concentration exists for which the effective diffusion coefficient 
reaches its maximum value. For a > 10, interactions monotonically reduce 
the effective diffusion coefficient as c increases. As in the previous case, the 
non-monotonic behavior manifested by the correction factor agrees well with 
the non-monotonic behavior of the effective diffusion coefficient observed in 
simulations and experiments^^—. 



3. THE EFFECT OF SURFACE DIFFUSION ON THE 
ADSORPTION-DIFFUSION INTERPLAY 

The case in which the adsorbed particles may diffuse on the surface of 
the sorbent can be analyzed by adding a diffusion term to Eq. ( Fl2|) .-. The 



evolution equation for c 2 takes the form 

F> r r) 1 

(23) 



^c 2 = fc a ciZ(c 2 ) - k d c 2 + 



D 2 (c 2 )-^-c 2 
ox 



where D 2 (c 2 ) is the diffusion coefficient at the surface. 

An estimation of the contribution of surface diffusion to the adsorption- 
diffusion interplay can be performed by assuming that D 2 ~ L 2 /td where td 
is the characteristic diffusion time and L a characteristic length of the sorbent 
surface. In first approximation, the diffusion term in ( 123]) can be expressed as 



d_ 

dx 



D 2 (c 2 )-^-c 2 



~ — . (24) 



This type of approximation is frequently and succesfully used in transport 
processes problems to examine and estimate, for instance, viscous or thermal 
effects near to a surface. 

At low coverage, when the adsorbed phase is dilute, one can make the 
appproximation Z(c 2 ) ~ 1 — £c 2 . Eq. (123]) reduces to 

^c 2 = k a c 1 (l-£c2)-k e J f c 2 , (25) 

where k e J — kd (l — Tf) 1 ^ 1 ). This equation leads to the adsorption isotherm 

c 2 = 0( Cl ) = , (26) 

where now a = k e J /k a < a. Therefore, in the dilute case the influence of 
surface diffusion on adsorption is to reduce the desorption rate. As a conse- 
quence, the values of the effective diffusion coefficient of the system D e ff(c) 
become shifted to lower concentrations. 
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In the nondilute case, both Z(c 2 ) and f 2 must include non-linear contribu- 
tions. The expression of Z(c 2 ) is given in Eq. (I2"0j) . According to Eq. (JSJ), 
the activity coefficient can be written in terms of the non-ideal contribution 
of the chemical potential: f 2 = exp(Afi n i/RT). A virial expansion of A/i n j 
at the lower order in concentration leads to the relation f 2 — exp((3c 2 ), where 
(3 = B 2 j RT has dimensions of volume and B 2 (T) is the second virial coefficient 
which must take into account excluded volume effects and direct interactions. 



Using the previous expressions for Z(c 2 ) and f 2 (c 2 ), one obtains the gener- 
alized Langmuir equation 

JU 2 = k aCl (1 - ic 2 + bc 2 2 ) - kfc 2 , (27) 

where now the effective desorption rate is k e J = kd [l — T d 1 ^d, 1 (1 + fi c 2)\ ■ 
The coefficient /3 can be positive or negative depending on the temperature. 
In the non-equilibrium stationary case one obtains from Eq. (I27j) 

a - 7 + £ci 



c = Ci + 



1 ± Vl - 4(ci - J)(6 Cl + /3 7 )/(« - 7 + ^ci) 2 • (28) 



/?7 + 2&C! 
where 7 = r^ l k~ l . 

In Figure HI we show the nonideal correction factor A s d(c) for £ = 1 in 
the case of equilibrium adsorption with surface diffusion in terms of the total 
concentration c for low, Fig. HJi), and high temperatures, Fig. 2]m). As in Fig. 
2, interactions induce a non-monotonic behavior but in the present case, for 
a < 1 and low temperatures (/3 < 0), they lead to an enhancement of diffusion 
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FIG. 4. The nonideal correction factor for the case of surface diffusion A S( ^(c) in 
terms of the total concentration c for b = (1 + a) 2 /4, J = 0, 7 = r^/c" 1 = 0.25 and 
different values of a = kd/k a and (3. InUJi) the solid lines represent the behavior 
of A s d(c) for equilibrium adsorption whereas dashed lines correspond to A n j(c) for 
{3 = —1.5cm 3 (low temperatures) and a = 0.25, 0.85, 1.45, 3.45 and 13.45 (from 
red to purple). The nonmonotonic behavior of the effective diffusion coefficient is 
more pronounced for a < 1. ii) shows A s d for (3 = 0.5cm~ 3 (high temperatures) 
and a = 0.45, 1.5, 3.0, 7.5 and 15.0 (from red to purple) and £ = 1cm 3 . For the 
present model, high values of (3 > may induce a discontinuity of the correction 
factor. 



for large concentrations due to the attractive nature of direct interactions of 
the particles in channel. At low concentrations, diffusion is significatively 
lower than in the case of Fig. 2. When the concentration is higher and the 
temperature is such that (3 > 0, excluded volume effects cause suppression of 
diffusion, as observed in the figure. 
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4. CONCLUSIONS 

In this article, we have analyzed the combined effect of diffusion and adsorp- 
tion in porous materials by adopting an effective medium description based on 
non-equilibrium thermodynamics^. 

We have treated the case of very narrow pores, present in some porous 
materials, microfiuidic devices and in biological channels, in which particles 
move essentially along the direction of the pores and may occasionally be- 
come attached to their surface. The fact that the pores are narrow, makes 
it possible to adopt a quasi-onedimensional description in terms of an aver- 
age concentration along the channel and a kinetic mechanism governing the 
attachment-detachment of the particles to the surface of the channel. The 
model proposed considers the combined effect of these two mechanisms and 
analyzes particle transport by means of effective diffusion in which the dif- 
fusion coefficient depends on the adsorption kinetics described in turns by a 
generalized Langmuir equation. 

The analysis of the effective diffusion coefficient reveals the existence of a 
multiplicity of different regimes controlled by different adsorption-desorption 
kinetics and different coverages of particles at the surface. We have quantified 
the change of the diffusion coefficient and found a non-monotonic behavior 
of this quantity due to excluded volume interactions. These results are in 
qualitative agreement with previous ones obtained from numerical models. 

The model proposed could be extended to the case of different shapes of 
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the pores by considering entropic barriers emerging from the variation of the 
cross section of the pores. This generalization can be performed following 
the approach presented in Refs.— ^ The resulting model may represent a 
simple and useful tool to analyze particle transport in nanoporous materials 
and biological channels. 
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